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£_) ' In December 1907, Minkowski expressed the Maxwell equations in the very beautiful and compact 4- 

Q" 1 dimensional form: lor / = — s, lor F* = 0. Here 'lor', an abbreviation of Lorentz, represents the 4- 

I ■ dimensional differential operator. We study Minkowski's derivation and show how these equations general- 

ly r%' ize to their modern premetric form in the framework of tensor and exterior calculus. After mentioning some 

applications of premetric electrodynamics, we turn to Minkowski's discovery of the energy-momentum ten- 
sor of the electromagnetic field. We discuss how he arrived at it and how its premetric formulation looks 
like. 
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1 Introduction 

In September 1908, Minkowski [40] presented a 4-dimensional geometrical version of special relativity 
theory; for the historical context, one may compare Walter l63l . A bit earlier |39l , in December 1907, 
he gave a preliminary version, and he expressed Maxwell's theory of electrodynamics in a 4-dimensional 
form. For the Maxwell equations he found, as pointed out by Damour |7] in his historical analysis, the 
very compact representation 



lor/ = -s and lorF* = 0. 



(1) 



Here lor^ := d/dxh are the components of the 4-dimensional differential operator. The excitation / and 
the field strength F are represented in Maxwell's nomenclature b>o 
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with the imaginary unit i := \/— 1. The rows and columns of the matrices in (|2]l are numbered from 1 to 
4. The 4-dimensional electric current is denoted by s. The star * used in (Q])2 is a duality operator, that is, 



* E-mail: hehl@thp.uni-koeln.de 
Minkowski took D = e, H = m\ E = E, B = M, that is, for the excitation / ~ (e, m) and for the field strength F ~ (E,M). 
Even though this seems logical and appropriate (e, E correspond to electric and m, M to magnetic), the Maxwellian notation is so 
deeply ingrained in our subconsciousness that any competing notation seems to have no chance. 
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^hk := \thkimFim, with ihkim = ±1, as the totally antisymmetric Levi-Civita symbol and 61234 = +1; 
moreover, the summation convention is applied. Minkowski demonstrated ||39l that the equations in (|TJ 
are covariant under arbitrary Poincare (inhomogeneous Lorentz) transformations. 

For vacuum (in the "aether", as Minkowski said), excitation and field strength, in suitable units, are 
related by 

f = F or (f hk ) = (F hk ). (3) 

Accordingly, Minkowski found the correct 4-dimensional representation of the Maxwell equations in the 
context of special relativity theory. 

Later it became clear, perhaps for the first time to Einstein 1111 . that the Maxwell-Minkowski equations 
(HJ can be generalized such that they become covariant under general coordinate transformations (diffeo- 
morphisms) and metric independent. Then, only the generalization of the constitutive relation (01 remains 
metric dependent. 

In the calculus of exterior forms, the 4-dimensional Maxwell equations for the excitation H and the 
field strength T eventually reaqj 



dH = J and dT = . 



(4) 



Here d is the exterior differential. The Minkowskian scheme (Q]) foreshadows in its compactness the exte- 
rior calculus version in ©. This point was stressed by Damour (7). The equations in © are premetriqj 
and, accordingly, also valid in this form in general relativity and in general relativistic field theories with 
torsion and curvature [48 1. The 2-forms of the excitation H = | Hijdx 1 A dx 1 and of the field strength 
T = i Tijdx 1 A dx 3 , with i,j = 0, 1, 2, 3, decompose in 1+3 according to |[T6l 



H = -HAdt + V and T = E Adt + B . (5) 

The signs in (O are not conventional. They are rather dictated by the Lenz rule, see Itin et al. I25ll26l . If 
we put the latter two equations into matrix form in order to be able to compare them later with the results 
of Minkowski, we find 



H = (H tj ) = 
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(6) 



Here V a := \e abc V bc and B a := \e abc B bc , with the 3-dimensional Levi-Civita symbol e 123 = +1, that 
is, T> and B, as required by their definitions in physics, are 3-dimensional 2-forms, whereas E and H are 
1 -forms. Note that V a and B a are vector densities. We stress that the identifications in (O are premetric, 
they are valid on any (well-behaved) differential manifold that can be split locally into time and space. 

The excitation H is a 2-form with twist (sometimes called odd 2-form) and the field strength T a 2-form 
without twist (sometimes called even form 2-form). These are consequences of the definitions of H and 
T via charge conservation and the Lorentz force density, respectively. For this reason, in exterior calculus, 
the constitutive law for vacuum ("the spacetime relation") reads 

H = \ *T. (7) 



Here Ao = yeo/Vo ~ 1/377 SI is the vacuum admittance; in the units Minkowski chose, Ao = 1. The 
metric-dependent Hodge star operator * in (O maps a 2-form without twist into a form with twist, and vice 
versa. 



We take the notation of our book [161 . The only exception is that we denote the field strength by T in order not to mix it 
up with Minkowksi's F. For other presentations of Maxwell's equations in exterior calculus, see, e.g., Bamberg & Sternberg [5], 
Choquet-Bruhat et al. QQ, Delphenich (HOD, Frankel (H), Jancewicz EJJ, Lindell (34), Misner, Thome & Wheeler ED, Russer (5T1 . 
Thimng |60l . and Trautman 1611 . 

We call a structure premetric, if it can be defined "before" a metric is available (or a metric exists, but it is not used). 
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In this paper we want to trace the way from Minkowski's system (Q]i and (O with (01 to the premet- 
ric framework (01 and (O with the metric-dependent spacetime relation ^}. Subsequently, we analyze 
Minkowski's discovery of the energy-momentum tensor of the electromagnetic field and display its mod- 
ern premetric version. 

2 Minkowski's way to: lor f = — s, lor F* = 

Minkowski first introduced the fields / and F in Cartesian coordinates x,y, z and with an imaginary time 
coordinate it (he put c = 1), see (|2]). Afterwards he changed to the numbering of coordinates in the 
following way: x\ := x, X2 '■= y, X3 := Z, X4 := it. His metric is Euclidean ds 2 = dx\ + dx\ + 
dx\ + dx\ = ghkdxhdxk, with g^k — diag(l, 1, 1, 1). Sometimes he also seems to use the signature 
(—1,-1,— 1,-1). In any case, as long as Minkowski sticks to Cartesian coordinates — and this is what 
he was doing — there is no need to distinguish contravariant (upper) from covariant (lower) indices. He 
displays the conventional Maxwell equations in component form. The inhomogeneous Maxwell equation 
reads 



d/12 d/13 d/14 
8x2 dx3 8x4 



Sl, 



d/21 d/23 d/24 

q ~r t S2 , 

OX 1 OX3 OX4 

"5 — + "5 — + "5 — = s 3 , (8) 

OX 1 OX2 OX4 

dfn d/42 d/43 

a ' n + a — S4 : 

0x1 0x2 0x3 
and the homogeneous one 

8F 3 4 , 9F 42 , aF 23 



9x9 9x3 <9a; 



= 0. 



3 wa-4 



ap43 + ap^ aFk 

dxi 8x3 8x4 

dF 24 , dF 41 , dF 12 



dxi 8x2 8x4 

dF 32 , 9F 13 , 9F 21 



= 0, (9) 

= 0. 



dx\ 8x2 dx3 

Using modern notation, including the summation convention, see Schouten 1521 . we can rewrite (0 and 
©as 

dfhk , dFhk . dFki dFi h 

— — = s h and — h h — — = . (10) 

OX k OXi dx h ox k 

Here h, k, ... = 1, 2, 3, 4. Minkowski also introduced the dual of Fhk according to Fj£ k := ^ehklm^lrm 
with the Levi-Civita symbol ihkim = ±1, and 61234 = +1. Explicitly, we have 
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(11) 



V -B x -By -B z J 



Then we can put the homogeneous equation in a form reminiscent of the inhomogeneous one (apart from 
the sources): 

r)F* 

^P*=0. (12) 

ox k 
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Subsequently Minkowski develops a type of Cartesian tensor calculus with a 4-dimensional differential 
operator called 'lor' (abbreviation of Lorentz). He introduces ordinary (co)vectors (space-time vectors of 
the 1st kind), like Xh and lor/j := -^—, and antisymmetric 2nd rank tensors (space-time vectors of the 2nd 
kind), like fhk and Fhk- Then, he can represent (TTOl ii and (fT2l symbolically as 

lor/ = -s and lorF* = 0, (13) 

respectively, see ([TJi. Using his Cartesian tensor calculus, Minkowski showed that equations ( fT3l are 
covariant under Poincare transformations. 

Einstein and Laub lfl3l thought that Minkowski's presentation of electrodynamics was very demanding 
from the point of view of the mathematics involved. Therefore, they essentially rederived Minkowski's 
results in a purely 3-dimensional formalism. 

3 Premetric electrodynamics in tensor calculus 

The next step in the development of the Maxwell equations occurred immediately after Einstein's funda- 
mental 1915 paper on general relativity and even before his big survey paper on general relativity would 
appear. Now Einstein was in command of tensor calculus in arbitrary coordinate systems. Einstein ifTTl 
observed that Maxwell's equations in vacuum can be put in a generally covariant form by picking suitable 
field variables. 

Einstein used the metric ds 2 = g^dx^dx" with signature (1, — 1, — 1, —1), here fi, v, ... ~ 0, 1, 2, 3. 
Einstein IfTTl wrote Maxwell's equations a£] 

dF p dF a dF Tp _ _ r— nmnVfi „ &$» 



= 0, r v = y/=ggr>tr p F afi , -£— = j». (15) 



dx T dx p dx a ' dx 

The field strength F pa is a tensor, the excitation J*"* a tensor density. In his "Meaning of Relativity" lfl2l . 
in the part on general relativity, he even picked the letter </> for the field strength (instead of F) apparently 
in order to stress the different nature of (j> and ^. 

The Maxwell equations ( TTBT li and ([T5l h are generally covariant and metric independent. Since in general 
relativity the metric g is recognized as the gravitational potential, it is quite fitting that the fundamental 
field equations of electromagnetism do not contain the gravitational potential. The equations in (TT~5T > are 
specifically valid in Minkowskian spacetime. Then, of course, the metric occurring in ([TBI l? is flat, i.e., 
its curvature vanishes. Hence by transforming Minkowski's lor / = — s and lor F* = to curvilinear 
coordinates, one can equally well arrive at (TT3T >. 

The gravitational potential only enters equation ([T5I19. We call this equation the spacetime relation — 
it is the "constitutive law" of the vacuum. Needless to say that having understood that there exists a 
way to formulate Maxwell's equations in a generally covariant and metric-independent manner, going 
back to ( [Tol l would appear to be an anachronism. Einstein's argument in favor of his formalism, which 
is similar to the Minkowski version ( [Tol l, and not following ( TTOb i and (fT2l . was that with his formalism 
one can derive the energy-momentum tensor of the electromagnetic field in a much more transparent way. 
This is certainly true. Still, it was Minkowski who discovered the energy-momentum expression l39l by 
performing products between the two space-time vectors of the 2nd kind / and F. 



Einstein used subscripts for denoting the coordinates x, i.e., x T etc. Moreover, we dropped twice the summation symbols S. 
Einstein's identifications, which were only worked out for vacuum, read 
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(14) 



Einstein discussed already earlier the general covariant form of Maxwell's equations, see |10| . However, the presentation was a bit 
different from (T5J and also by far not as transparent as the one in 1111 . Probably for this reason Einstein republished these ideas. 
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But the real decisive step was left to Kottler (29]. He derived the Maxwell equations from the con- 
servation laws of electric charge and magnetic flux. The charges in a 3-dimensional volume and the flux 
lines piercing a 2-dimensional area are conserved. These two independent conservation laws are related to 
counting procedures and thus are independent of any length or time standard. It is for this reason that the 
two Maxwell equations (fT3T >i and dl~5l >^ emerge as a priori independent structures and free of a metric. It 
is then left to the constitutive law to link the excitation ^ u with the field strength F a p. 

In dT3b we displayed the Maxwellian system in the way Einstein wrote it. Changing now to modern no- 
tation, see Schouten lf52l and Post l47l , and using the conventions of lfT6l , the premetric Maxwell equations 
reac0(i,j,... = O,l,2,3) 



djH* 3 
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(16) 



The antisymmetrization bracket is defined according to [ijk] := ^(ijk—jik+jki — [-•■•) i.e., a plus (minus) 
sign occurs before even (odd) permutations of the indices ijk. Here H^ are components of a tensor density, 
whereas jFy are those of a tensor. For Tij we have the identification ©2 and for H^ := ie y Hm, with 
Hki identified according to ©1. Thus, explicitly^ 
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-CijklJ > then ( [ToT l can be displayed in the more symmetric 
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diiTjk] = . 



(19) 



Today it is not only of academic interest to put Maxwell's equations in a general covariant form. On 
the surface of a neutron star, for example, where we might have strong magnetic fields of some 10 10 tesla 
and a huge curvature of spacetime — we are, after all, not too far outside the Schwarzschild radius of the 
neutron star of some 3 km — the Maxwell equations keep their form ( fTol ). The metric of spacetime enters 
the vacuum relation, see (fT5]>9. There is no longer place for only a Poincare covariant formulation of 
Maxwell's equations a la (1101) , A less extreme case is the Global Positioning System. Nevertheless, also 
the GPS rules out Poincare covariant electrodynamics, see Ashby JH. 

3 . 1 Constitutive relation 

The system of the 4 + 4 Maxwell equations (TToT l for the 6 + 6 independent components of the electromag- 
netic field H lJ and Tki is evidently underdetermined. To complete this system, a constitutive relation of 
the form 



#« = #«(^h) 



(20) 



has to be assumed. The constitutive relation ( |20| | is independent of the Maxwell equations and its form can 
be determined by using experimental results. For general magnetoelectric media, including the vacuum as 



Post [47| denotes H*i by © A " / . In Obukhov's investigations on the energy-momentum tensor 1421 , the H*i is simply named 
H'i , that is, the " is dropped. 

If we abandon Minkowski's imaginary time coordinate and go over to x 1 , with i = 0, 1, 2, 3, then excitation and field strength 
in Minkwowski's framework read 
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special case, we assume with Tamm l57l[36l . v. Laue 11321 . and Post l47l the /ocajj and linear premetric 
relation 

H ij = ±x iJM fki, (21) 

where x is & constitutive tensor density of rank 4 and weight +1, with the dimension [\] = 1/resistance. 
Since both H'i and Tki are antisymmetric in their indices, we have x 1 ^ — — X = — X • An 
antisymmetric pair of indices corresponds, in four dimensions, to six independent components. Thus, the 
constitutive tensor can be considered as a 6 x 6 matrix with 36 independent components. 

One can take the premetric Maxwell equations ( fT6l ) together with the premetric ansatz d2"TT i as an elec- 
trodynamic framework ("local and linear premetric electrodynamics"), see also Matagne 1381 . Therein one 
can study the propagation of electromagnetic waves in the geometrical optics approximation. It turns out, 
see J50l[T6l . that the premetric, totally symmetric Tamm-Rubilar tensoiQ 

Q m (X) := ^ e mnpq e rstu X mnr(i X 3lps]k X l)qtu (22) 

— a new structural element in theoretical physics — controls the light propagation via a generalized Fresnel 
equation Q % ^ kl qiq^q^qi = 0. The q L are the components of the wave vector of the lightjj If birefringence is 
forbidden, then the light cone emerges |30l , that is, the metric g^ up to a (dilation) factor |fl9l . 

Let us come back to x l3kl ■ As a 6 x 6 matrix, it can be decomposed in its tracefree symmetric part (20 
independent components), its antisymmetric part (15 components), and its trace (1 component). On the 
level of x i m is decomposition is reflected in 

ijkl _ (1) ijkl + (2) ijkl _j_ (3) ijkl _ ^23) 

36 = 20 © 15 © 1. 

The third part, the axion part, is totally antisymmetric and as such proportional to the Levi-Civita symbol, 
(3)^ijki ._ ^[ijki] _ a e ijki j\j ote that a is a pseudoscalar since e l: > kL has weight — 1. Therefore, the weight 
of e yfe ' is essential information. The second part, the skewon part, is defined according to t 2 )^*^ := 
i^ijkl _ -^khjy jf the constitutive equation can be derived from a Lagrangian, which is the case as long 
as only reversible processes are considered, then ( 2 ^ l i kl = 0. The principal part ^^ kl fulfills the 
symmetries W^ kl = t 1 )^-^ and ( 1 )^[*-' fc '] = 0. The premetric constitutive relation now reads 

H ij = i ( {1) x ijkl + {2) X m + ol e llkl \ Th ■ (24) 

If we assume the existence of a metric gij, then in vacuum the constitutive tensor reduces to 

ijkl \ / / ik jl il jk\ /oc\ 

X =MV-9{9 9 ' -9 9 ' >■ (25) 

This formula is, up to a factor, determined uniquely by the symmetries of the principal part. On substitution 
of (f25t into ( ETl i. we find the vacuum relation 
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= \oV=gg* k g Jl F k i, or ■- ^jmnV^gg mk g nl F M , (26) 

{X) v ijkl 



compare (TT~5b-j. Note that (f25b belongs to the principal part ^\ 



A nonlocal generalization of this law is being investigated by Mashhoon, see his review 1371 . 

Parentheses around indices denote total symmetrization, that is, {ijkl) := -^{ijkl + jikl + jkil +■■■}. The covariant 
epsilon-system is marked by a hat: e,jkl = ±1,0; eoi23 = +1- 

This scheme was generalized to Abelian gauge and string theory by Schuller, Wohlfahrt, et al., see the articles |53, 49] and the 
references given there. Also in these generalized schemes a Tamm-Rubilar type tensor plays a decisive role. 
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In order to compare $7M with experiments, we have to split it into time and space parts. As shown 
in [[TBI ITTl in detail, we can parameterize the principal part by the 6 permittivities e ab = e ba , the 6 
permeabilities p, ab = p ba , and the 8 magnetoelectric pieces j a b (the trace vanishes, j c c — 0) and the 
skewon part by the 3 permittivities n a , the 3 permeabilities m a , and the 9 magnetoelectric pieces s a b . 
Then, the premetric constitutive relation (f24b can be rewritten as 



V a = (e ab - e abc n c ) E b + ( 7 % + s b a - 8 a b s c c ) B b + aB\ 

H a = (p-i - e abc m c ) B b + (- 7 b Q + s a b - 8 b a s c c ) E b - a E a . (27) 

Here e abc = e abc = ±1,0 are the 3-dimensional Levi-Civita symbols, with e 123 = £123 = +1. As 
can be seen from our derivation, a is a 4-dimensional pseudo (or axial) scalar, whereas s c c is only a 3- 
dimensional scalar. The cross-term j a b is related to the Fresnel-Fizeau effects. The skewon contributions 
m c , n c are responsible for electric and magnetic Faraday effects, respectively, whereas the skewon terms 
s a b describe optical activity. Equivalent constitutive relations were formulated by Serdyukov et al. 11541 . 
p. 86, and studied in quite some detail. 

3.2 Some applications of premetric structures 

The theory of a possible skewon part of ( l24b lies in its infancy, see IF2T1I441 . The axion part was widely 
neglected in the literature, if not even put to zero right away l47l . As we can see from (124-b and d271 l. the 
axion piece obeys 

<«>*<-!<,«««*, o, \ll_fA. (28) 



• Tellegen 1581 1591 introduced the gyrator, with a = const, as a new device in the theory of two- 
port networks that "rotates" the field strength (B, E) into the excitation (T>, Ti) — that is, "voltages" into 
"currents" — and vice versa, see ifTTl . 



• In the axion-electrodynamics of Wilczek 1651 . the axion piece as a field is added to 

H lJ = (a„ V~ 9 g lk g 31 + \a{x)e^ kl \ F kl , (29) 



see also Itin 1241 . Moreover, kinetic terms for the axion are added in order to make the axion a propagating 
field. 



• Lindell & Sihvola [35. 561 . again for a = const, defined the perfect electromagnetic conductor (PEMC) 
via (l28l . Hopefully the PEMC can be realized by suitable metamaterials, see Sihvola 1531 . Lindell & 
Sihvola arrived at the PEMC when studying electrodynamics in its 4-dimensional version. Then the minus 
sign in H = —aE emerges naturally. Note, T> = —sB and H = —sE, with s := s c c , see (1271) . represent 
a (spatially) isotropic skewon that is only a 3-dimensional scalar. Thus, the relative sign in the second 
column of (1281 is decisive. 



• On the basis of experiments in solid state physics, we were able to show [20] that in Cr2C<3 a phase can 
occur (below the Curie temperature) that obeys the relation (written in terms of exterior forms) ^H = a.T, 
with a k, 10~ 4 Ao. This proves that H is a form with twist and T one without twist thereby demonstrating 
the consistency of the system dH = J, dJ- = and ^H = aJ-. 

4 Premetric electrodynamics in exterior calculus 

One of the principal goals of exterior calculus is to put equations, in the case under consideration, Maxwell's 
equations and the constitutive law, into a manifestly coordinate independent form. In a way, in exterior 
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calculus, the principle of coordinate covariance is built-in. If we handle the 2-form H, we handle a geo- 
metrical quantity which does not depend on coordinates at all. Only if we are interested in its components 
Hij, with H = ^ Hijdx 1 A dxK then these components are, of course, coordinate dependent. Thus, the 
Maxwell equations as such, namely dH = J and dT = are just relations between forms, coordinates do 
not play a role. Naturally, such a calculus expresses the universality of Maxwell's equations in a particular 
impressive way. 

Let us derive the Maxwell equations in modern language ab ovo in order to liberate the physics of 
Maxwell's equations from historical ballast, see IfTSI for details. Electrodynamics is based on two conser- 
vation laws. First, we have electric charge conservation (first axiom). The charge-current density 3-form 
J is closed dJ = and, by de Rham's theorem, also exact, 

dH = J, (30) 



with the excitation 2-form H as "charge-current potential" [62]. This is already the inhomogeneous 
Maxwell equation. 

The current J is a 3-form with twist "since an electric charge has no screw-sense," cf. Schouten l52l . As 
a consequence, the excitation H also carries twist or, in the words of Perlick [46| (my translation), "...one 
must understand the excitation as a form with twist if one wants that the charge contained in a volume 
always has the same sign, independent of the orientation chosen." 

We decompose the 4-dimensional excitation H into two pieces: one along the 1 -dimensional time t and 
another one embedded in 3-dimensional space (a, b — 1,2, 3). We find (see lfl6l ) 

H = -H A dt + V = -H a dx a Adt+ -V ab dx a A dx b (31) 

which, if put in matrix form, yields the identification ©l . Substitution of OTb into ( f30b leads to the 
3-dimensional inhomogeneous Maxwell equations, 

dV = p, dU-V = j, (32) 

i.e., to the Coulomb-Gauss and the Oersted-Ampere-Maxwell laws, respectively. The underline denotes 
the 3-dimensional exterior differential and the dot differentiation with respect to time. 

With charge conservation alone, we arrived at the inhomogeneous Maxwell equations ( 13 Ot . Now we 
need some more input for deriving the homogeneous Maxwell equations. The force on a charge density is 
encoded into the axiom of the Lorentz force density (second axiom) 

f a = (e a \T)AJ. (33) 

Here e a is an arbitrary (or anholonomic) frame or tetrad, a basis of the tangent space, with a = 0, 1, 2, 3, 
and J denotes the interior product (contraction). In tensor language, we can express (l33l as /< = Tik J k . 
The axiom (l33l should be read as an operational procedure for defining the electromagnetic field strength 
2-form T in terms of the force density f a , known from mechanics, and the current density J, known from 
charge conservation. According to J33I . the field strength T turns out to be a 2-form without twist. Its 1+3 
decomposition reads 

T = ENdo + B = E a dx a f\dt + -B ab dx a A dx b , (34) 

which yields the identification 4(6)2- 

Magnetic flux conservation is our third axiom. The field strength T is closed, 

dF=0. (35) 

This is the homogeneous Maxwell equation. Split into 1+3, we find 

dE + B = 0, dB = 0, (36) 

Copyright line will be provided by the publisher 



adp header will be provided by the publisher 



that is, Faraday's induction law and the sourcelessness of B. The laws ( I36l l are also premetric. The 
Lenz rule, as the reason for the relative sign difference between the time derivatives in (|32] >9 and (l36lli. is 
discussed in Il25ll26l . 

We have then the Maxwellian se^j 

dH = J and dT = . (37) 

If decomposed into components, we recover the tensor calculus formulas in (TT~9b . This closes our con- 
siderations on the appropriate form of Maxwell's equations in 4-dimensional spacetime. For vacuum, the 
Hodge star maps the 2-form T without twist into the 2-form H with twist according to 

H = A *T (38) 



— this is equivalent to (126l l. The Hodge star, as applied to 2-forms in 4 dimensions, depends only on the 
conformally invariant part of the metric, see, e.g., Frankel lfl4l and Kastrup 11281 . The constitutive relation 
for axion-electrodynamics simply reads 

H = X Q *T + aT, (39) 

that is, for the PEMC we have simply H = ouF. For matter we have the relation 

H = k{T) (40) 

as analog of ( f20b . If the operator k is local and linear, we find for its components 

1 

2' 



K .. kl — -2.. ^rankl ^1) 

•Hj — c 2j ran A V* 1 -/ 



5 Minkowski's discovery of the energy-momentum tensor 

Minkowski's "greatest discovery was that at any point in the electromagnetic field in vacuo there exists a 
tensor of rank 2 of outstanding physical importance. . . . each component of the tensor E p q has a physical 
interpretation, which in every case had been discovered many years before Minkowski showed that these 16 
components constitute a tensor of rank 2. The tensor E p q is called the energy tensor of the electromagnetic 
field." 

These words of Whittaker J64l . Vol. 2, pp. 66 and 67, stress the central importance of this discovery. 
Einstein wouldn't have been able to derive the field equations of gravity of general relativity without being 
aware of this important quantity that figures as source of his (and Hilbert's) field equation of gravity. 
Einstein's reason for picking the scheme ( fT5l l. and thus (fTST l with (l26lii. was that this allowed him to derive 
the energy-momentum tensor of the electromagnetic field in a more transparent way than Minkowski. 

Minkowski starts with his "field vectors" / and F, as specified in (0, and multiplies them as matrices 
according to fF. The product matrix he decomposes, without giving any motivation, into a tracefree piece 
S (that is, with tr S = 0), which he calls a spacetime matrix of the 2nd kind, and a trace piece C 1, with the 
unit matrix 1 and with C := —jtr(fF), which Minkowski identifies correctly as the Lagrangian density 
i (H B — D-E) of the electromagnetic field. He wrote, if we formulate it compactly in the notation of 
Cartesian tensor calculus, 

Shk = IhiFik - -ShkfimFmi , that is, S hh = . (42) 



Damour (7) discussed Maxwell's equations in exterior calculus in his footnote 1, He found 8f = s and 5 * F = 0. The 
codifferential is defined according to <5 := *d*. For a p-form ijj we have **i/> = (— l) p ~ 1 »/>. For the homogeneous equation we find 
*d**F = — * dF = or dF = 0, which coincides with J371 ? . Thus F can be identified with the field strength T. We apply the 
star to Damour's inhomogeneous Maxwell equation and find * (<5/) = **d*f = d * f = S := *s or d * f = S. This equation 
only agrees with (37t i , provided we put H = */. Consequently, in Damour's exterior calculus, / loses its operational interpretation 
as electromagnetic excitation. Of course, the codifferential <5 is not required for the formulation of the Maxwell equations (37). The 
exterior differential d is sufficient. 

Copyright line will be provided by the publisher 



10 



Friedrich W. Hehl: Maxwell's equations in Minkowski's world 



In order to display the reality relations appropriately, he wrote his spacetime matrix as 

/ Sn S12 S13 

_ / q \ _ S21 S22 S23 
" — \i->hk) — c c c 

<J31 ^32 -J33 

y 5*41 S42 543 

With the help of ( f42b , Minkowski identifies explicitly the Maxwell stress (momentum flux density) 



Si 4 \ 




( x x 


Y, 


Z. r 


-*T X \ 


S24 




Xy 


Yy 


Zy 


-iTy 


S34 




x z 


Y z 


z z 


-iT z 


S44 / 




\ -iX t 


-iY t 


—iZ t 


T t J 



(43) 



X x = - (D X E X — D y Ey — D Z E Z 

Xy = D y E x + H x B y , 

Y x = D x Ey+H y B x etc., 

the Poynting vector (energy flux density) 



H X B X — HyB y — H Z B Z ) , 



T„ 



H,E,, 



H y E z 



etc. 



and the electromagnetic energy density 



T t = \ (D X E X 



DyEy 



D Z E Z + H X B X 



HyBy 



H Z B Z ) 



(44) 



(45) 



(46) 



He did not name the electromagnetic momentum density X t , Y t , Z t — which, by Lebedev 11331 . had already 
been discovered experimentally in 1901 . Note that Minkowski's energy-momentum tensor is asymmetric in 
general and thus, because of its tracelessness, has 15 independent components (for the algebraic properties 
of energy-momentum currents, see Itin J23l ). 

We leave this historical track and turn first to tensor calculus: Starting with the Lorentz force density 
fi = TijJ k , substituting the inhomogeneous Maxwell equation d!6) i. and using (TToTk we arrive at the 
Minkowski energy-momentum tensor densit>[|J 



7? = -8]T kl H kl -T lk W k 



% k 



0. 



(47) 



together with /, = djTj + Xt and the auxiliary force density Xj = j [{diTjk) H^ k — ^FjkdiH^]. 

i S k T 2 , with T 2 := % m T m \ a formula that 



Incidentally, already Minkowski proved |39| that %° Tf 
was repeatedly rediscovered. 

Vox vacuum, H % i = Xq \J—gg lk g : ' 1 Fkh and thus, in Cartesian coordinates [gij = diag(l, — 1, —1, —1)], 
the auxiliary force density Xi = 0. Then % k := gkj%? = Iki is symmetric and conserved, see |fl6l for 
details — and we recover the result of Einstein [TT). At the same time, under those two conditions (vacuum 
and Cartesian coordinates), Minkowski's formula K = lor S re-emerges, where K is the 4-dimensional 
Lorentz force density and S Minkowski's energy-momentum matrix. 

In exterior calculus we can apply an analogous procedure and find the energy-momentum 3-form as 



1 



[T A{e a \H)-H A{e a \T)] , #* A S Q = , 



(48) 



where e a denotes the frame (tetrad) and d a the coframe, with e a \ {)P = <5g. The expression J48] ii translates 
into d47i ii via E Q = T a fi *-dp. Accordingly, the formulas d47b and (f48b incorporate Minkowski's energy- 
momentum of the electromagnetic field in a modern notation. 

Apparently Minkowski's energy-momentum tensor is defined for matter, in particular also for moving 
matter. Abraham HI 12 introduced a symmetric tensor instead, since the symmetry was believed to be 



11 If we use HV = \e^ M H kl , then 7? = \ e> klm (H ih F lm - F ik H h 
inspired by Minkowski's formalism, displayed the energy-momentum tensor as T : 



It is remarkable that Ishiwara 1221 in 1913, 

-i{[[H*F*]]-[[FH]]}. 
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important from the point of view of angular momentum conservation, amongst other reasons. The dispute 
between the supporters of the Minkowski and the Abraham tensor is lingering on till today. Obukhov l42l 
reviewed this problem, partly following 11431 [161 l45l . Obukhov 11421 based his discussion on a general 
Lagrange-Noether framework. He recovers the Minkowski tensor as the canonical energy-momentum and 
shows that the balance equations of energy-momentum and angular momentum are always satisfied for an 
open electromagnetic system because of the asymmetry of the canonical tensor. Moreover, Obukhov is able 
to formulate an "abrahamization" prescription for symmetrizing a given asymmetric energy-momentum 
tensor, provided a timelike vector is available (a velocity of a medium, for example). Obukhov considers 
an ideal fluid with isotropic electric and magnetic properties interacting with the electromagnetic field. 
Thereby Obukhov resolves the perennial Abraham-Minkowski controversy and assigns the appropriate 
places for the Minkowski as well as for the Abraham tensor. 



6 Discussion 

Compare, in vacuum, Minkowski's system 

{lor/ =-s, lorF* =0; / = F} , (49) 

with the system in exterior calculus 

{dH = J, dT = 0; H = *T} . (50) 

The first system is Poincare covariant, the second one generally covariant. 

Equation (1491 is a special case of (l50l . Since in Minkowski's formalism lor F* = in its version dTOb ? 
is already generally covariant, we rewrite it tentatively as (lor* ) F = and identify lor* — ► d. Then, the in- 
homogeneous equation becomes (lor*) f*=sordf* = s, and we arrive at {d f* = s, dF = 0; / = F}, 
or at ( TSUl l. provided we identify H — > /* and T — ► —F. Roughly we could say that in the transition from 
J49l to (T50b the star migrated from the Maxwell equations to the constitutive relation for the vacuum. Of 
course, the system d49l is now history. The set d50l l. which was derived from a reasonable axiomatics of 
Maxwell's theory, has superseded it. Still, we find it remarkable how powerful Minkowski's formalism is 
and how closely d50b resembles d49l . 

And in exterior calculus Minkowski's energy-momentum matrix Shk finds its particular compact ex- 
pression as energy-momentum 3-form 

Z a :=^[TA(e a \H)-HA(e a \T)] . (51) 

Summing up: The Maxwell equations dH = J, dT = 0, the Lorentz force density f a = (e Q J J-) A J, and 
Minkowski's energy momentum 'tensor' dSTb represent the framework of premetric electrodynamics. As 
such, they are a modern generally covariant expression of Minkowski's 4-dimensional Poincare covariant 
theory of electromagnetism. 
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